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ABSTRACT 

In general relativity static gaseous atmospheres may be in hydrostatic balance in the 
absence of a supporting stellar surface, provided that the luminosity is close to the 
Eddington value. We construct analytic models of optically thin, spherically sym¬ 
metric shells supported by the radiation pressure of a luminous central body in the 
Schwarzschild metric. Opacity is assumed to be dominated by Thomson scattering. 
The inner parts of the atmospheres, where the luminosity locally has supercritical 
values, are characterized by a density and pressure inversion. The atmospheres are 
convectively and Rayleigh-Taylor stable, and there is no outflow of gas. 
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1 INTRODUCTION 


Several system s with super-Eddington luminos ity have been 
reported (e.g., iMcClintock fc Remillard 200^}jthe “LMC 
transient” A0535-668 i Bradt fc McClintock 19831 ) being oft 
discussed. Recently, interest in very luminous neutron stars 
has been revived with the discovery of a clear example of a 
neutron star with super-Eddington flux in the guise of the 
ultraluminous 1.37 s pulsar NuSTA R J095551-I-6940.8 in the 
nearby galaxy M82 llBachetti et al.ll2014l ) . It is clear that for 
some accreting neutron stars, radiation pressure may exceed 
the pull of gravity close to their surface (at least in a certain 
solid angle if the radiation is beamed). 

In this paper we report the presence of a new type of 
atmospheric solutions in general relativity (GR) for neutron 
stars radiating at nearly Eddington luminosities. These so¬ 
lutions are qualitatively different from the ones obtained 
in Newtonian physics. In the classical solutions the atmo¬ 
spheric density increases monotonically as the stellar surface 
is approached, and this remains tru e in previously obtained 
atmo spheric solutions in GR (e.g., IPaczynski fc AndersonI 
Il986l) . where the atmospheres in question have always been 
supported at their base by the stellar surface. Here, we re¬ 
port solutions in which the atmosphere has the form of a 


shell suspended above the stellar surface. The maximum of 
pressure is attained on a surface enclosing a volume in which 
the star is contained (together with some space around it), 
and the atmospheric density and pressure drop precipitously 
on both sides of that surface, i.e., in the radial direction away 
from the star (as is usually the case), but also in the radial 
direction towards the star. 

This paper discusses the simplest case of optically thin, 
spherically symmetric atmospheres in the Schwarzschild 
metric, which admit of analytic solutions. Optically thick at¬ 
mospheres require numerical treatmen t of radiative transfe r, 
and these will be reported elsewhere (IWielgus et akllioish . 

We consider atmospheres consisting of pure ionized hy¬ 
drogen. All the results are given for the Schwarzschild space- 
time. The luminosity of the star will be parametrized by the 
ratio of the luminosity observed at infinity to the Eddington 
luminosity. 


— Too/TEdd, 


( 1 ) 
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with the standard expression for the latter, Z/Edd = 
A-KGMrripc/aT. The stellar radius will be denoted by R*, 
it is taken to correspond to a canonical neutron star, e.g., 
R* « 5GM/c^, with M « I.AMq being the stellar mass. 
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2 BALANCE BETWEEN GRAVITY AND 
RADIATION PRESSURE 


In the Schwarzschild metric, the stellar luminosity (which we 
take to be constant in time), as measured by a local observer 
at radius r, declines with the radial coordinate distance r as 


L{r) = 


rc^ ) 


( 2 ) 


and a static balance between gravity and radiation pressure 
owing to Thomson scattering can only be achieved at one 
radius, at which 

L{r) = LBdd(l - Rs/r)-'/", (3) 

where Rs = 2GM/rc? is the Schwarzschild radius (e.g., 
IPhinne^^l 19871 : [Stahl et al.ll2013h . That radius can be written 
as 


r = Rs/ (l — A^) = Regs- (4) 

The sphere at Regs has been called the Eddington cap¬ 
ture sphere (ECS) bv IStahl et alj (120121 '). These results have 
been derived rigorously within the mathematical formalism 
of Einstein’s general relativity, and they can be intuitively 
understood in terms of the different dependence on redshift 
of the luminosity and of the “gravitational” acceleration. 
Eq. dU reflects the dependence of luminosity at infinity on 
two redshift factors (1 -I- z)~^, one for the energy of the 
photons, and one for the rate of their arrival. However, the 
acceleration of a static observer in the Schwarzschild metric 
scales with only one factor of {1 + z), where the redshift is 
defined by 

l + z(r) = (l-Rs/r)-i/^ (5) 


Hence, the condition for local balance between accelera¬ 
tion and the momentum flux of photons (multiplied by the 
Thomson opacity) is given by Eq. (|3]), and contains one fac¬ 
tor of {1 -\- z), which introduces into the balance condition 
a dependence on the radial coordinate. This is in contrast 
with Newtonian gravity (which can formally be recovered in 
the limit Rs —^ 0, or z —>■ 0), where the balance condition 
has no radial dependence, and therefore is either satisfied, 
or not satisfied, equally at all distances from the source. 

Numerous authors have shown that test particles ini¬ 
tially orbiting the star (at various radii) may settle on the 
spherical surface at r = Regs > R», provided that 

(1-Rs/R.)^^" < A < 1, (6) 


their angular momentum having been removed by radiation 
drag jBi ni et al.ll2009l: ISo k Oh et ^l201ll : IStahl et '^120121 : 
iMishra &: Kluzniak 2014 1. In fact, any point on the ECS 
is a position of stable e quilibrium in the radial direction 
dAbramowicz et^ll99d l. and neu tral equilibrium in direc¬ 
tions tangent to the ECS surface dStahl et al.l 12012 1. Note 
that the redshift z(r) attains the following value on the ECS, 



Figure 1. Density profiles for the isothermal solutions with 
T = 10^ K. Each curve corresponds to a different value of the 
luminosity parameter: A = 0.8, 0.85, 0.9, 0.95, 0.97, 0.99 from left 
to right (as labeled at the maxima). 


the fluid is optically thin. This will allow us to decouple 
the temperature of the fluid from the radiation field and 
to find analytic solutions of the atmospheric structure. It 
is clear that the condition of hydrostatic equilibrium will be 
achieved if the pressure gradient is in the radial direction and 
its magnitude compensates the imbalance between radiation 
pressure and gravity. In particular, the pressure must reach 
a maximum value at r = Regs, and must decrease in both 
directions, towards and away from the star, as the distance 
to the ECS grows and so does with it the imbalance between 
radiation pressure and gravity. 

For atmospheric temperatures of no more than a few 
keV, we can safely assume that the energy density of the 
fluid is given by its baryon rest mass energy density pc^. In 
other words, one can neglect the contribution of pressure p 
and the internal energy e to the energy density, pc^ ^ p + e. 
The equation of hydrostatic equilibrium then becomes 


1 dp 
p dr 


GM 

r^ (1 — Rs/r) 


1 - A 



( 8 ) 


where the last term reflects the radiation pressure owing 
to Thomson scattering, with the luminosity of the star 
scaled by the Eddington luminosity per Eq. O- In the limit 
Rs/r —S' 0 one recovers the classical Eddington balance at 
A = 1, while in general the pressure gradient vanishes for 
r = Regs, cf. Eq. (H. 

Eq.® is readily solved for a polytropic or an isothermal 
atmosphere, as it can be integrated over the redshift to yield 


2EGS = z{Recs) = 1/A — 1. (7) 

We now take the step of replacing test particles on the 
ECS with a fluid. For simplicity, let us assum^B that 

^ The optically thic k case will be reated in a separate paper 
llWielgus et abllioiih . 


[ = ln(l + z) — Xz. (9) 

J pc 

For an ideal gas equation of state, p = phsT/{prtip) with ks 
the Boltzmann constant, rrip the proton mass, and 
the mean molecular weight, we obtain the following analytic 
solutions. 
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3 ISOTHERMAL ATMOSPHERE 


For an isothermal atmosphere at temperature T, the integral 
in Eq. m yields 

In f ^ = ln(l + z) — \z. (10) 

with z defined in Eq. ((5]). In the GR case, the maximum 
value of density is attained at the ECS, assuming a lu¬ 
minosity in the range given by Eq. ®, and the density 
can be taken to be proportional to this maximum value 

pmax — P(zeCs), 

P(Z) = Pm.. X [A(l + 2)](M-P==)/(fcBT) 

X exp{[l — A(1 -I- z)](pmpC^)/(kBT)j . (11) 


Thus, as expected, the density decreases in both directions 
away from the ECS, giving the atmosphere a shell-like form. 
The same result, given as a function of the radius, is 

p(r) = Pm.. X [A(l - Rg/r)-i/2](M-pc=)/(feBT) (12) 

X exp |[1 - A(1 - Rs/r)~^^^](pmpC^)/(kBT)j . 


However, far from the star this solution is unphysical, as the 
mass in an isothermal atmosphere would be infinite—note 
that at 2 = 0, i.e., as r —>■ oo, the density tends to a finite 
value, p^poo= Pm.. X [Aexp(l - ^ 

The Newtonian Eddington limit is recovered with a dif¬ 
ferent normalization from the one in Eqs. (HU-dm. To first 
order in z, from any of the three Eqs. ([Tni)-([T^ one gets, 
p(r) = Poo X exp [pmpGM(l — A)/(kBTr)]. The classical Ed¬ 
dington limit, dp/dr = 0, or p = const for an isothermal 
atmosphere, is now obtained for A = 1, while the Newto¬ 
nian limit of an exponentially decaying atmosphere holds 
for |r — ro| < ro 


p(r 


po exp 


pnip GM{1 — A) 

/crT rn 


{r - ro) 


(13) 


with arbitrary constants po and ro of dimension density and 
length, respectively. 

Although we have assumed Thomson scattering, this is 
not strictly necessary. In fact, all the above considerations 
in Section[3l as well as Eq. remain valid if we replace the 
Thomson value, ctt, of the photon electron scattering cross 
section by a more general cross-section o{T), e.g., the Klein- 
Nishina formula, provided that we make the substitution 


A —>■ Ai = Aai, with ai = <j{T)/gt- 

We stress once again that unlike the decaying atmo¬ 
sphere of the Newtonian limit in Eq. GSl), the GR solution 
in Eqs. ([T0ll-([T^ is not monotonically decreasing with ra¬ 
dius everywhere, but instead is monotonically increasing for 
r < Recs up to the maximum on the ECS, and is only 
monotonically decreasing for r > Regs, as illustrated in 
Fig. (O. 


4 POLYTROPIC ATMOSPHERE 

For a polytrope, p = Kp^, the specific enthalpy w = 
(p/p)r/(r — 1) following from Eq. dS]) or Eq. dS]) is 

w{z) = wbcs + <? {In [A(l -I- z)] — A(1 -I- z) -|- 1} (14) 

and it attains its maximum value, wecs = ui(2bcs), at r = 
Regs, see Eq. ©• 



Figure 2. Temperature profiles for the pol 3 Ftropic solutions. Each 
curve corresponds to a different value of the luminosity parameter: 
A = 0.8, 0.85, 0.9, 0.95, 0.97, 0.99 from left to right. Note that 
the temperature difference on the y axis is in units of 10® K. 


Making use of the ideal gas equation of state we obtain 
a corresponding solution for the temperature 


T(r) — Tm.. = 
F — 1 prUpC^ A 
~ ^) 


(15) 




-In 




also attaining its maximum value on the ECS. Note that for 
any given value of the luminosity. A, the temperature falls 
of on both sides of the ECS in a universal way (Fig. [2} inde¬ 
pendently of the value of the maximum temperature Tm.., 
which is just an additive constant in Eq. m- Of course 
the temperature cannot go to negative values, so the plot 
must be cut off at T — Tm.. = —Tm..- Fig. [2] corresponds to 
Tjnax = 10® K, i.e., a hot ~ 100 keV corona. 

The density scales with its arbitrary value at the ECS, 
pm.. = p(Regs), which is proportional to the total mass of 
the atmosphere; 

p{r) = po |^ln(l - Rs/r)“^^® - A(1 - Rs/r)“^/® -|- Aj , 

(16) 

with p(Regs) = Po[A — In A — 1]^'^^^ 

This being a polytrope, the atmosphere has a sharp 
edge, the enthalpy, temperature and density vanishing at 
the top of the atmosphere. In fact, when the condition of 
Eq. di) is satisfied, there are two “tops” of the atmosphere 
(Fig. [S]), one at r+ > Regs, with z+ = z{r+) < zecs, and 
one at r_ < Regs, with 2 _ = 2 (r_) > zbgs- The redshift 
at the top of the atmosphere is given by the condition 

{In [A(l -I- 2 ±)] - A(1 -I- z±) -I- 1} = -wbgs/c®. (17) 


For maximum temperatures not much larger than a few keV, 
the atmospheric profile is quite symmetric with respect to 
the sign of r — Regs (Fig. [3]), and the height of the atmo¬ 
sphere is proportional to the speed of sound at its “base,” 
Cs = y(r — 1 )ibegs- 

Indeed, for H = \r± — Regs| -C Regs, 

H y2A® c.c 
Regs y/V — 1 ’ 

where Uk = \/ GM/Rbgs is the Keplerian orbital velocity 
at the ECS. This is reminiscent of the result for accretion 
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Figure 3. Density profiles of polytropic Thomson-scattering at¬ 
mospheres for A = 0.9. For a luminosity this large, the atmosphere 
is well separated from the neutron star surface (at R^/M 5). 
Temperatures from outside in (more extended atmospheres to less 
extended atmospheres) are Tmax = 5 ■ 10^ K, 5 ■ 10® K, 5 ■ 10® K, 
5 • lO'^ K, respectively, with = 1/2 and F = 5/3 in all cases. The 
density maxima are at -Recs = 10.5GM/c^. 



Figure 4. Same as Fig.[3l but for A = 0.95. The density maxima 
are at Recs = 20.5GM/c^. 


disks, where Hjr^CsjvK (|Shakura &; Sunvaevll 19731 ). how¬ 
ever the ECS atmospheres are more extended, H here being 
enhanced by a factor of cjvK’ In terms of the temperature 
at the base of the atmospheric shell, this reads 

1/2 / U rj^ \ 1/2 

firUpC^ J 


H 


2y 


Recs l-AMr-l 


2r 


For ksTinax = 2 keV, and Loo ~ L^dd the shell thickness 
works out to be JL rs 2-10~^(1 — For a neutron star 

this gives H « 10(1 — A) “^m, i.e., H = 1 km for A = 0.9, and 
H = 4 km for A = 0.95. One can see that as the luminosity 
increases, the thickness of the levitating atmospheric shell 
increases more quickly than the radius of the ECS. 


5 DISCUSSION 

We have shown that luminous compact stars (e.g., neutron 
stars) in Einstein’s general relativity may have atmospheres 
which are detached from the surface. This may have applica¬ 
tions to astrophysical phenomena of those compact objects 
in which the luminosity is close to its Eddington value, ei¬ 
ther on a quasi-permanent basis or in a transient state. The 


former case includes highly luminous neutron stars accret¬ 
ing matter from low-mass stellar companions, such as the Z- 
sources (which include the brightest X-ray source in the sky, 
Sco X-1). One would expect quasi-spherical coronae in such 
sources to be described by our solutions with implications 
for the hard X-ray spectra in these sources, and possibly also 
for their quasiperiodic modulations. The latter include X-ray 
bursts, some of which exhibit the phenomen on of radius ex¬ 
pansion ie.g.. I^rohmaver fc BildstenI 120061 1. However, here 
we have only considered optically thin atmospheric shells, 
whereas the atmospheres of X-ray bursters are optically 
thick. For this reason we are postponing a detailed discussion 
of bursters to an other paper, where o ptically thick solutions 
will be reported (I Wielgus et al.l[2015li . 

The atmospheres presented here can be thought of as 
an (extreme) example of density and pressure inversion. At¬ 
mospheric density inversion owing to variations of opacity 
has been discussed previously in various contexts. In stellar 
structure theory discu ssion of densi t y inve rsion goes back to 
IChandrasekhai lll936h . see also lSerJ 1 194l| j. In contemporary 
papers such layers are often invok ed in the context of wind 
losfl and stellar pulsations (e.g., iTownsend fc MacDonaldl 
| 200 ^ . 


The possibility of gas pressure inversion was pointed 
out by Erika Bohm-Vitense, and explained as a result of a 
change of sign in the effective gravitational acceleration, geft, 
with ffeff being the difference between the true acceleration 
of gravity and the acceleration owing to radiative forces, 
which are pro portional to the pro duct of opacity and the 
radiative flux (IBohm- Vitensel[i958h . This is indeed the cause 
of the pressure inversion in the atmospheres presented in 
our paper. However, there are several important differences 
between the results obtained previously and those reported 
in this work. 

First, the inversion layers familiar from stellar structure 
theory are caused by variations of opacity, while our results 
are of a different origin—the (Thomson) opacity in our so¬ 
lutions is uniform, and the density (and pressure) inversion 
arises solely owing to effects of general relativity. More gen¬ 
erally, its presence is related to purely geometrical effects 
that allow the radiative flux to have a different functional 
dependence on the distance from the source than the accel¬ 
eration of gravity. 

Second, we are presenting optically thin atmospheres, 
whereas the stellar discussion occurred for optical depths 
greater than unity. 

Third , in the optically thick case considered by 
IJoss et ^ (Il973l i density inversion occurs in regions where 
the luminosity is near-Eddington (although still subcritical), 
the pressure is dominated by radiation, and the temperature 
gradient is superadiabatic. While the atmospheres discussed 
in our paper do require near-Eddington luminosities, none of 
the other conditions is met. Density inversion occurs only in 
the region of supercritical luminosity, the radiation pressure 
is zero, as the radiation freely streams through the optically 


^ In LBV and hypergiant stars “the most striking propert[y] is 
the strong density inversion in the outer layers, where a thin 
gaseous layer floats upon a radiatively supported zone. This is 
due to the peak in the opac ity which for ces supra-Eddington lu¬ 
minosities in some layers.” llChios il ll99^ 
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thin gas, and the temperature gradient is subadiabatic for 
the isothermal atmosphere, of course, as well as for poly¬ 
tropic atmospheres with F < 5/3. Gas-pressure inversion 
occurs at supercritical luminosities. The region of pressure 
inversion coincides with that of density inversion in the op¬ 
tically thin atmosphere s prese nted here. 

Fourth, IJoss et ^ (Il973ll . in their discussion of opti¬ 
cally thick inversion layers in hydrostatic eq uilibrium, re¬ 
ferring to the work of iLucv fc SolomonI lll970li observe that 
their own “results do not alter the well known conclusion 
that mass-outflow from the surface must result if L > L/rit, 
where L and Tcrit a-re the photospheric values” of luminosity 
and its critical value. We note that our atmospheric solu¬ 
tions do “alter the well known conclusion,” as the levitating 
atmospheric shell is in hydrostatic equilibrium even though 
L > Z/Edd in (a part of) the optically thin shell, i.e., the lumi¬ 
nosity is supercritical above the photosphe re of the underly¬ 
ing lu minous star. This is possible because ILucv fc SolomonI 
lll97nll used Newtonian gravity while our solutions are valid 
in GR. 

We conclude with a brief discussion of stability. The at¬ 
mospheres presented in this paper are convectively stable, 
as shown by the appropriate Schwarzschild criterion, in the 
isothermal case and for polytropes with F < 5/3. This result 
can be intuitively understood by realizing that by the law 
of Archimedes a hot parcel of gas is initially accelerated in 
the direction opposite to the gradient of pressure, hence the 
pressure-inversion layer is as stable as the layer of radially 
decreasing pressure. A poly tropic atmosphere with F = 5/3 
is marginally stable to convection. The atmospheres are also 
Rayleigh-Taylor stable as the denser layers are always “be¬ 
low” (in the sense of the direction of effective gravity) the 
less dense layers. Investigation of the normal modes of the 
atmospheres will b e presented elsewhere. Her e, we only note 
that (as shown by lAbarca fc KluzniakI [20151) the radiation 
drag efficiently damps oscillatory motion in a class of fun¬ 
damental modes, just a s it does the rad ial and azimuthal 
motion of test particles llStahl et al.ll201^ . 
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